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General Solution to Gradient Induced Transverse and Longitudinal Relaxation of 

Spins Undergoing Restricted Diffusion 
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We develop an approach, by calculating the autocorrelation function of spins, to derive the mag- 
netic field gradient induced transverse {T2) relaxation of spins undergoing restricted diffusion. This 
approach is an extension to the method adopted by McGregor. McGregor's approach solves the 
problem only in the fast diffusion limit; however, our approach yields a single analytical solution 
suitable in all diffusion regimes, including the intermediate regime. This establishes a direct con- 
nection between the well-known Torrey's slow diffusion result and the fast diffusion result. We 
also perform free induction decay measurements on spin-exchange optically polarized ^He gas with 
different diffusion constants. The transverse relaxation profiles are compared with the theory and 
satisfactory agreement has been found throughout all diffusion regimes, fn addition to the transverse 
relaxation, this approach is also applicable to solving the longitudinal relaxation (Ti) regardless of 
the diffusion limits. It turns out that the longitudinal relaxation in the slow diffusion limit differs 
by a factor of two, compared with that in the fast diffusion limit. 

PACS numbers: 33.25.-(-k 32.60.-(-i 34.10.-l-x 51.20.-|-d 
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I. INTRODUCTION 

Longitudinal relaxation (Ti relaxation) and transverse 
relaxation (T2 relaxation) are the most important param- 
eters in Nuclear Magnetic Resonance (NMR), Magnetic 
Resonance Spectroscopy and Magnetic Resonance Imag- 
ing (MRI). Many factors can contribute to Ti and T2 
relaxations. It is well known that, for liquid or gaseous 
samples, such as spin polarized '^He, diffusion in a non- 
uniform magnetic field can cause both Ti and T2 relax- 
ations. In this manuscript, 1/rp and I/T2 represent 
the relaxation rates solely due to diffusion in a magnetic 
field gradient. The transverse relaxation due to diffusion 
in free space was first solved by Torrey [l|. He general- 
ized the Bloch equation by adding a diffusion term and 
unveiled that the envelope of the transverse component 
of the magnetization decays as 



A{t) 



eM-^Dl^Gh''), 



(1) 



where D is the diffusion constant, 7 is the gyromag- 
netic ratio of the spin and G is a constant gradient. In 
practice, most diffusion happens in confined spaces, and 
the stochastic diffusion process is restricted, which makes 
the problem more complicated. Robertson [2] solved the 
Bloch- Torrey equation [l| in restricted geometries, by im- 
posing boundary conditions to the equation. His approx- 
imate analytical solution showed that the envelop decays 
exponentially with a constant relaxation rate. This re- 
sult is valid when ADt >> L^, where L is the distance 
between the boundaries. This limit is also known as fast 
diffusion limit or motional averaging regime, where spins 
have moved across the geometry many times in a time 
period t and therefore any fluctuation in the magnetic 
field averages out and a faster diffusion actually reduces 
the relaxation. Neuman [^ solved the same problem by 
calculating the accumulated phases of spins with the as- 



sumption that the relative phase distribution of spins is 
Gaussian in both the slow diffusion and fast diffusion 
limits. His slow diffusion result reproduces the free dif- 
fusion result, Eq. ([I}, and the fast diffusion result is the 
same as that of Robertson. In the intermediate regime, 
the Gaussian Phase Approximation (GPA) fails. How- 
ever, it is crucial to quantitatively understand the inter- 
mediate regime because many experiments have shown 
edge enhancement phenomena in the slow diffusion and 
intermediate regime, which was recognized later as a lo- 
calization regime [j|. People have observed that, when 
water diffuses in microscopic structures, the MRI signal 
is enhanced at the edge of the structure [^-Q ■ Saam et al. 
have also showed a similar edge enhancement effect, using 
hyper pol arized '^He gas in cells with dimensions of about 
1 cm [i3|. This effect is ascribed to the more restricted 
diffusion at the boundary, which lessens the relaxation, 
and was first described quantitatively by De Swiet |ll| . 
using Airy functions. Airy functions have been shown 
to be the eigenfunction of the Bloch- Torrey equation in 
the intermediate regime [T3- Axelrod also showed that 
although GPA fails in the intermediate regime, it can be 
used to interpolate the result in this regime, which turns 
out to be close to the exact solution [lj|. More detailed 
discussion on the restricted diffusion in various limits can 
be found in a review article JJJ and references therein. 

Despite of the widely used GPA method. Gates et al. 
used the second order time-dependent perturbation the- 
ory and carried out an expansion of spin density ma- 
trix to obtain both longitudinal and transverse relax- 
ation rates for a spherical cell '15'|. Their results works 
only in the fast diffusion regime, and they further di- 
vided the fast diffusion regime into two limits: the high 
pressure limit ujoR'^/SttD >> 1 and the low pressure 
limit ujqR'^ /SttD << 1, where luq is the Larmor preces- 
sion frequency and R is the radius of the spherical cell. 
These two limits can be thought of as the characteristic 



spin precession time ti = 2tt/ujo being much shorter or 
longer than the characteristic diffusion time in the cell 
Trf = R^ /AD, respectively. McGregor [1^ used Redfield 
theory discussed in \\J\ to solve the same problem also 
in the fast diffusion regime. Redfield theory is a general- 
ized treatment of the second-order time-dependent per- 
turbation theory. It establishes a set of differential equa- 
tions obeyed by the spin density matrix. Therefore, it is 
closely related to the treatment of Gates et al. [l^l . By 
calculating autocorrelation function of spin under differ- 
ent geometries, McGregor was able to obtain transverse 
relaxation rates in the fast diffusion limit for different 
geometries, including slabs, cylinders and spheres. 

In this manuscript, we make an extension to McGre- 
gor's approach, which yields an analytical solution to 
the transverse magnetization suitable for all diffusion 
regimes. When ADt « L^, this solution reproduces 
Torrey's free diffusion result; when time ADt « i^, it is 
in the intermediate regime. In these two regimes, edge 
enhancement effect is also observed. Eventually, when 
ADt >> L^, it is in the motional averaging regime and 
a peak located at the center of the frequency spectrum 
is observed. We also performed Free Induction Decay 
(FID) measurements on polarized '^He gas to verify the 
theoretical results. By changing the number density of 
the gas, the observed transverse relaxation happens in 
different diffusion regimes. When the decay envelopes are 
compared to the theoretical predictions, they are found 
to be in good agreements, especially in the intermedi- 
ate regime. In addition to the transverse relaxation, our 
approach can also be used to calculate longitudinal relax- 
ation in different regimes. We found that the longitudinal 
relaxation rate l/Tp in the slow diffusion limit is twice as 
fast as that in the fast diffusion limit. As diffusion in the 
fast diffusion regime is more restricted, it could explain 
this factor of two difference. In this manuscript, we solve 
the problem in ID for clarification purpose. However, 
it can be easily extended to 3D with complex geome- 
tries since one only needs to calculate the corresponding 
probability density function. Once the density function 
is known, the relaxation rate can be calculated readily, 
which makes this approach suitable for numerical simu- 
lations of complex geometries. 



II. REDFIELD THEORY FOR MAGNETIC 
FIELD GRADIENT-INDUCED RELAXATIONS 

For simplicity, let spins diffuse in a cubic cell with 
length L. A non-uniform magnetic field is applied along 
the z direction. At time t — 0, we track a spin starting 
at x'^. As time evolves, the expected position of the spin 
will change due to the diffusion process. Therefore, we 
use {x'{t)) to represent the expected position of spin at 
some later time t and has the property that (a;'(0)) = x'^. 
Since the field is non-uniform over the box, the spins also 
see fluctuating magnetic fields during diffusion. The fiuc- 
tuating field B' can be treated as a perturbation to the 



zeroth order mean field B{{x')) by taking Taylor expan- 
sion around x = {x'), 

B' ^ B{x{t)) Biix'it))) ^ VB ■ {x{t) - {x\t))). (2) 

As described in [13] , Redfield theory gives solutions to the 
problem with fluctuating magnetic fields. In our case, 
the applied field is constant in time; however the time 
dependence appears because of the diffusion process. The 
equation of motion for the transverse and longitudinal 
components of the spin can be written as |17| 

I {St) ^ I ((5.) + ^ {Sy)) = 7((St) x B)- 
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Jb'{uj) + Jb'H +Jb,(0)W5t) (3) 



dt 



{S, 



-1 



Jb'M + Jb'Au:) (S,) 



(4) 



where St is the transverse component of spins, B = 
Bz{{x'))z, ijj — jBz{{x')), the Larmor precession fre- 
quency of spins at field strength Bz{{x')), and Jb' is 
defined as 111 



Jb'^ (w) 
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(5) 



The bar denotes an ensemble average of the autocorre- 
lation of the perturbed magnetic field. Jb', (0) has the 
similar definition with B'^ replaced by B'^ and a; = 0. 
The first term in Eq. ([3]) describes the precession of the 
spin under the field Bz{{x')), and the second term gives 
the transverse relaxation rate l/Tj^; whereas Eq. ([4]) 
describes the longitudinal relaxation. 

It should be noted that Eq. ([5]) was originally written 
in [ITJ as an integral from to infinity. It was argued 
that the autocorrelation of the magnetic field would van- 
ish quickly after a critical time Tc, and consequently, in- 
tegration from to infinity introduces negligible errors as 
long as i >> Tc- For the case considered here, Tc can be 
defined as t^ = {L/2)'^/AD. When the diffusion is slow, 
the above approximation is invalid, and Eq. ([5]) must be 
used. By utilizing Eq. ([2]) and realizing the fact that dif- 
fusions in X, y and z directions are independent of each 
other, Jb' (ijj) can be re-written as 



,dB' 



,dB' 



,dB' 



(6) 
where Jx{<-^) is expressed as 



Ulo) = / {x{t - r) - {x'{t))){x{t) - {x'{t)))e-'^-dT, 
Jo 

(7) 
and similar for Jy{uj) and Jz{^)- In general, Jx{^) is 
time-dependent and the relaxation rates are also time- 
dependent. However, it will be shown that, in the fast 



diffusion limit, Jx{oj) approaches some constant values 
and Jxi^j) = Jy{uj) — Jz{uj). Hence, one can define the 
longitudinal relaxation rate 1/Tp as 



Tf 



-fWVB,\' + \VBy\')^[Uu)] 



(8) 



where /9(a;,i|xo,io) is the conditional probability density 
function of spin at position x at time i, if the spin is at 
xq at time to [to = t — t). It is apparently that Eq. pH)) 
satisfies the condition (a;'(0)) ~ x'^. 



where SR [•] indicates taking the real part. Hence, in the 
fast diffusion limit, the transverse relaxation rate 1/TJp 
can also be defined as 



1 



1 



T§ 2Tf 



-7'ivB,r j,(o). 



(9) 



A. Magnetic Field Gradient-induced Transverse 
Relaxation Rate IjT^ 



The expected position of spin {x'{t)) at time t can be 
computed as 



L/2 

{x'{t))= j xp{x,t\x'^,0)dx, 

-L/2 



In order to compute Jx{0), one needs to calculate po- 
sition autocorrelation function, which, by definition, can 
(10) be expressed as 

I 



rL/2 .L/2 

{x{t - t) - {x')){x{t) - (x')) = / {xo~ {x)){x~ {x'))p{x,t\xa,to)p{xo,to\(x') ,0)dxdxo. (11) 

J-L/2J-L/2 
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Due to the symmetry of the cubic cell, p{x,t\xo,to) can 
be found by decomposing the 3D diffusion equation into 
ID problem, subject to the boundary condition 

d 

— p(x,t|xo,io)U=±L/2 = 0, 

and the initial condition 

p{x,to\xo,to) = d{x- xo). 
The solution is found to be 



p{x,t\xQ,to) = 



1 2 
L + L 

L ^ 

n=2,4 



E_ .i^t^ct m:x mrxo 

e ~?^sm(— — )sm(— — ) 



n=l,3 

„2_2 , 



mrx , , mrxo . 
cos(— — )cos(— — ). 



(12) 



By substituting Eq. ([T2|) into Eq. (fTTj) . one obtains 
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{x{t -r)- {x')){x{t) - {x')) = — ^ - exp( 



r! = l,3 
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, n-K Dt. , knDto. ,kiT {x 
exp( r^ — ) exp( — ) cos( 
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fcTT 



(13) 



Instead of substituting Eq. P^ into J^ directly, one can 
first take the slow diffusion or fast diffusion limits on it. 
It can be shown (see Appendix) that, in the slow diffusion 
limit 4Di << (i/2)2, Eq. ^ is simplified to 



{x{t - r) - {x')){x{t) - {x')) w 2D{t - t). (14) 

Ja;(0) is then calculated as 



J^(0) = / 2D{t - T)dT = Dt^ 
Jo 



(15) 



This is also true for Jy{0) and Jz(0). As Jx,y,z{0) is time- 
dependent, the relaxation rate is also time-dependent. 
One has to substitute Jx,y,z{0) into Eq. (jH]) directly to 
compute St^ 

I {St) = 7((St) xB((x')))-(^+7'|Vi?.pi?t2) (St) . 

(16) 
If one assumes that the gradient in one direction 
dBz/dx = G dominates, then |Vi3z| ~ G and B^ — 
Bz{0) +G{x') w B,{0) + Gx'^. {x') « 4 because spins 
are localized in the slow diffusion limit and the expected 
position at time t will not change significantly from their 
initial positions. Integrating upon time, Eq. (|16l) be- 
comes 



{St) = 5oexp[z7(B40)+Gx^)t]exp(- 



-i'G'Dt' 



t 



2rp' 



(17) 

This coincides with ID diffusion result derived by Torrey 
[l|. Torrey's derivation assumes free diffusion without 
boundaries, which is equivalent to the slow diffusion in a 
confined volume, as boundaries are not present to spins 
when diffusion is slow. 

In the fast diffusion hmit, ADt » (£/2)^, all the ex- 
ponential terms containing t in Eq. (I13p vanish, so does 
{x') term, which can be proved easily by taking the limit 
ADt/L^ -^ oo in Eq. (ITU|) . Therefore, the only surviving 
term is the first one. 
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{x{t-T)-{x')){x{t)~{x'))^ J2 ^^e 



n=l,3 



Jx{0) in this case becomes 

/' y. "■" 
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Jx{0) 
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"dr 
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1201?' 



(18) 



(19) 



where the fact that ADt >> (i/2)^ is used again after 
the integration. Jy(0) and Jz(0) are the same as Jx{^) 
in this limit. Substituting them back into Eq. ([3]), one 
obtains 



{St) = 5*0 exp[ia;oi - (-— + 
111 
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12QD 



-)t]. (20) 



As {x'{t)) approaches zero in the fast diffusion limit, aU 
the spins precess at the same central frequency uoq — 



•yBziO), no matter where the spin is initially. This is 
also known as motional averaging regime. If gradient in 
one direction dominates the other two (IVB^I « G), the 
transverse component decays with a constant relaxation 
rate l/T"!^ given by 



1 



1 72L4G2 



2Tf 



+ 



120D 



(21) 



This result is also derived by McGregor [T^ and Robert- 
son [25 using GPA method. However, it will be shown 
in the discussion section that, by numerically calculating 
Jx (0) using Eq. ()13p without any approximation, one can 
obtain the frequency spectrum of the precession signal in 
various limits, which cannot be obtained from McGre- 
gor and Robertson's methods. As a constant gradient is 
applied, the frequency spectrum is actually a frequency 
encoded ID image. In the slow diffusion limit, peaks 
are observed at the edge of the geometry, known as edge 
enhancement; whereas in the fast diffusion limit, a reso- 
nance peak is observed at the center of the spectrum. 

In the intermediate region ADt « {L/2)'^, one has to 
substitute Eq. (|T3l) directly into Jx{0) and then calcu- 
late Eq. ([3]). Unfortunately, no concise analytical form 
of St{{x')) can be obtained in this regime. Hence, a 
numerical calculation of St{{x')) is performed and com- 
pared with Free Induction Decay (FID) measurements on 
gaseous '^He cells, and a good agreement in the interme- 
diate regime is shown, see Sec. HI. 



B. Magnetic Field Gradient-induced Longitudinal 
Relaxation Rate l/Tf 

In order to compute Eq. Q, one needs to obtain Jx{^) 
first. In the slow diffusion limit {ADt « {L/2)^), Eq. 
^fT^ should be used, 

Jx{io)= f 2D{t~T)e~"^''dT 
Jo 

^ 2D{1 - e-^^' - lut) 

Since, in most cases, uj >> 1, 5R(1 — e^*"^*) = 1 — cosujt 
is a fast oscillating function, which averages to 1. Con- 
sequently, substituting ^[Jx{uj)] into Eq. dS]), the longi- 
tudinal relaxation has an averaged decay rate as 

iVBxi' + lVByf 



1 2n 

— =^\\S/Bx\' + \S/By\')-^^2D 



Bl 



(23) 



where lo = -iB,{{x')) « 7^,(0) = 7S0. 

In the fast diffusion limit {ADt » {L/2)^), Eq. (H 
should be used to calculate Jx{tij), 



Jx{uj) 
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(24) 



Hence, the real part of Jx{^) is written as 

nU^)\ = ^^ E ^ ^,.n^...r. - (25) 



r! = l,3 



^L^' 



In the fast diffusion limit with high pressure {t^Jti - 
ujL'^ jZli^D >> I), the sum in Eq. ([25]) is simphfied to 



E 

n=l,3 



1 1 



n2a;2L4 gco^L'^-^ 



(26) 



and l/T^ in this hmit becomes 



(27) 



It is interesting to see that there is a factor of 2 difference 
between the slow diffusion 1/Tp and the high pressure 
fast diffusion l/T^ . More discussion on this topic is pre- 
sented in Sec. IV. 

In the fast diffusion with low pressure (rd/r; — 
L^L"^ JH-kD « I), the sum in Eq. ^^ becomes 



E 

n=l,3 



1 1 



n6 'K^D'^ 960152 ' 



(28) 



and the resultant ^/T^ in the low pressure limit is 

1 72^4 



Tf 1201) 



|VS,|' + |VBJ'). 



(29) 



This result is an analogy to the low pressure 1/rp de- 
rived by Gates et al. for a spherical cell geometry |15j . 
For geometries other than sphere and box, one only needs 
to recalculate Eq. p^ and the corresponding l/T^ can 
be obtained readily through steps illustrated above. This 
also applies to the transverse relaxation rate I/T2 when 
other geometries are considered. 



C. Magnetic Field Gradient-induced Resonance 
Frequency Shift 

In Eq. ([3]) , the imaginary part of the complex function 
Jb' and Jb' gives rise to the shift of precession frequency 

a; y 

diO, 



5uj 



-7 



\VB,\^ + \VB,/ 



3[Jx(c^)] 



(30) 



where 5 [•] means taking the imaginary part. In the slow 
diffusion limit, substituting Eq. ([22|) into Eq. (|30l) yields 



(31) 






It is interesting to note that, in the slow diffusion limit, 
the frequency shift increases linearly as a function of 
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Helmholtz coils 



FIG. 1: The apparatus for measuring the transverse relax- 
ation of ^He in the cylindrical measurement cell. "^He in the 
detachable cell has been polarized by SEOP before it is trans- 
ferred to the measurement cell. 



time, different from the t^ dependence in the transverse 
relaxation rate. In addition, Eq. (1311) does not depend on 
L as expected because the slow diffusion limit is equiva- 
lent to the free diffusion, in which spins do not see bound- 
aries. 

In the fast diffusion limit with high pressures, substi- 
tuting Eq. ([M| into Eq. ([50)) and taking the correspond- 
ing limit yields 



,,2 r2 



5lo 



^(|VB.|V|Vi?,|^) 



(32) 



In the fast diffusion limit with low pressures, it yields 

^ 17^07 B ^i^g |2 [yo |2) (33) 

20160£)2 ^' '■^\ \ v\ ) \ I 

These two results are analogies to the frequency shifts 
derived in [15| for a spherical cell. 



III. EXPERIMENTS AND RESULTS 

FID measurements have been performed to measure 
transverse relaxation of polarized '^He gas at 34.5 kHz. 
The '^He gas is polarized in a 2 inch diameter spheri- 
cal detachable cell via Spin-Exchange Optical Pumping 
(SEOP) technique. The cell is routinely filled with 760 
torr ■^He and 100 torr N2. After ^He is polarized, the de- 
tachable cell is connected to a measurement cell sitting 
at the center of the Helmholtz coils through a 2 mm in- 
ner diameter glass tube. The measurement cell is made 
of bare pyrex glass in a cylindrical shape with 2 inch for 
both the diameter and the length. Three intermediate 
volumes made of pyrex and a turbo pumping line are 
connected to the glass transfer tube between the detach- 
able cell and the measurement cell. The schematics of 
the experimental apparatus is shown in Fig. [TJ Before 
the measurement, polarized '^He atoms are allowed to dif- 
fuse into intermediate volumes first so that the number 
density of '^He in the detachable cell is diluted. The inter- 
mediate volumes are then isolated and the valve between 



the detachable cell and the measurement cell is opened. 
Consequently, the rest of '^He in the detachable cell can 
diffuse freely into the measurement cell. When the pres- 
sure in the two cells reaches equilibrium, the valve is 
closed and FID measurements are carried out. 

A pair of 6 inch diameter pickup coils in Helmholtz coil 
configuration is used, so that the coil has a rather uniform 
sensitivity over the entire cylindrical cell. Each coil has 
2000 turns of 30 AWG wires to maximize the signal. The 
precession signal at 34.5 kHz is lock-in detected and the 
envelop of the signal is extracted. The RF coil is a pair 
of saddle coils with a length of 3.5 inch, a diameter of 3 
inch and the opening angle is 120 degrees. Each coil has 
10 turns of AWG 22 wires. The axis of the pickup coil, 
the RF field direction and the magnetic holding field are 
perpendicular to each other. A RF pulse with a tipping 
angle ^ 20 degrees is sent to the measurement cell. A pair 
of gradient coils in Maxwell coil settings [lOl] is also added 
to provide a uniform field gradient of 2.3 mG/cm in the 
holding field or z direction. The background gradients 
are measured to be much smaller than this value and 
therefore ignored in the calculation shown below. 

In Fig. ^, we show the transverse relaxation mea- 
surements of '^He with pressures of 327 torr, 2.94 torr 
and 1.59 torr. The corresponding diffusion constants are 
3.5, 388 and 717 cm^/s, respectively (scaled from 1.5 
cm^/sec for 760 torr "^He, which is measured using Carr- 
Purcell-Meiboom-Gill method [20] )■ When using the ra- 
tio {L/2)'^/ADt to characterize the diffusion regime, it 
is clear that 3.5 cm^/s is in the slow diffusion regime; 
whereas 717 cm^/s is in the fast diffusion regime and 388 
cm^/s is in the intermediate regime close to fast diffu- 
sion limit. When diffusion is slow, the decay of signal 
is mainly due to the dephasing of the spin precessing at 
different frequencies. This can be understood by taking 
D as zero, so that the exponential term in Eq. p7|) is 
identical to one and a rough estimate of the overall signal 
Sail is simply integrating 5*^ over the entire cell, 



Sail 



-L/2 



L/2 



g.(^o+7G:.')rf2.' o, e*'^''*sinc(^^) (34) 



where sinc(x) is defined as sin(a;)/a; and the first zero of 
sine function is at ^GLt = 2tt. Using the real values of 
G and L, one can determine t = 26 ms. This value is the 
same as the experimentally measured signal vanishing 
time i « 26 ms for the case D=3.5 cm^/s, see Fig. ^. 

When diffusion is fast, spins will more or less precess at 
the same frequency (motional averaging) and the diffu- 
sion term dedicates the signal decay. As no simple form 
of St can be obtained in the intermediate regime, we 
numerically calculate St{{x')) using the general form of 
autocorrelation function, Eq. (IT^ . To evaluate Eq. ([T^ , 
we compute one hundred terms in each of the first three 
sums and four hundred terms in the last double sum. 
We also evaluate {x') up to 100 terms. Once St{{x')) is 
known, it is weighted by B{{x')), the sensitivity of the 
Helmholtz pickup coil at position {x'), and then inte- 
grated over the entire cylindrical cell to mimic the mea- 




- Signal 
-Simulation 



D=3.5 cm/s 



ijiikW^^ 



0.00 0.01 0.02 0.03 0.04 0.05 



C 300 

3 



CO 



-Signal 
-Simulation 






0.00 


0.O5 


0.10 


0.15 


500 












:,. 


ll 


Signal 

Simulation 




400 




- 






l-lll , 




D=717cm'/s 


300 


- 


rki 


1 II . 


- 


200 


(c) 


Hvi 


fHk* 


- 


100 


" 


' 




" 



0.05 0.10 

Time (s) 



FIG. 2: (Color online)Transverse relaxation measurements 
with different diffusion constants, (a) is in the slow diffu- 
sion regime and the sine Uke shape of the decay profile is due 
to the spin defocus, (b) is in the intermediate regime and (c) 
is in the fast diffusion regime. Simulation results are shown as 
dashed lines and compared to the measured decay envelops. 



sured FID signal. The simulation results, shown as red 
dashed curves, are compared with measured decay en- 
velops. The background noise of the FID measurement is 
around 150 arb. units. The simulation curves are shifted 
up by this amount to account for the background. Good 
agreements between measurements and simulations are 
found for the intermediate regime and the fast diffusion 
regime. One can also use Eq. (I^Tj) to predict T2 in Fig. 
(I2J;) as it is in the fast diffusion regime. The prediction 
yields T2 — 0.0589 s and an exponential fit of the data 
yields T2 — 0.0557 s, which is very close to the predic- 
tion. However, if Eq. ([?!]) is used to predict T2 in Fig. 
(I2}d), it overestimates T2 by 39% and the profile of the 
measured relaxation is somewhere between the sine and 
exponential. This shows that Eq. (I?T|) is inadequete to 
use in the intermediate regime, and one has to use the 
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FIG. 3: Frequency spectrum of St as a function of (a;') jL 
at f = 0.02 s. The length L of the cell is 1 cm. The solid 
line corresponds to D = 1000 cm^/s, which is in the motional 
averaging regime, and the dashed line corresponds to D = 1 
cm^/s, which is in the slow diffusion limit. 



non-approximated form of S'T((a;')) to do the calculation. 
In the slow diffusion regime, the numerical calculation 
correctly captures the time when signal vanishes, and it 
also exhibits a small bump at 0.04 s, due to the partial 
refocus of spin. However, this bump is not observed in 
the experiment. It is probably due to the fact that in 
the simulation, we only take into account the gradient in 
the longitudinal direction. In reality, although gradients 
in other directions are smaller than the longitudinal one, 
they still affect the precession frequency of each individ- 
ual spin. Consequently, spin refocus is disturbed and the 
small bump is smeared out. 



IV. DISCUSSION 

In Sec. Ill, we numerically calculate St{{x')) with- 
out any approximation. The time evolution of integrated 
St{{x')) is compared to the FID signal of polarized '^He 
gas and a good agreement is found, especially in the 
intermediate regime. It is also interesting to see how 
St changes as a function of position {x') when differ- 
ent diffusion regimes are considered. In Fig. [3J we show 
St as a function of {x'), which is also equivalent to a 
frequency spectrum due to the linear relationship be- 
tween uj and {x'), known as frequency encoding. The 
two curves shown in Fig. |3] are calculated at the time 
instant t = 0.02 s, with jG = 1000 rad/s-cm, L = 1 cm 
and D = 1 and 1000 cm^/s, respectively. The D = 1 
cm^/s case is in the slow diffusion limit, and two peaks 
close to the edges are observed. As diffusion is more re- 
stricted at the boundary, the diffusion induced relaxation 
is suppressed, compared with the relaxation at the cen- 
ter. In contrast, the D — 1000 cm^/s case is in the fast 
diffusion limit and only one peak centered at the mean 
frequency presents, which means most of spins precess at 
the same frequency and relax at the same rate, i.e. Eq. 
(I?!]) . These results show that the approach developed 



in this manuscript is able to capture all distinct behav- 
iors of the transverse magnetization in different diffusion 
regimes. 

In Sec. II B, it is shown that the longitudinal relax- 
ation rate l/T^p differs by a factor of 2 between the slow 
diffusion limit and the fast diffusion high pressure limit. 
A possible explanation is that when t is small, i.e. in the 
slow diffusion limit {ADt << (L/2)^), most of the spins 
do not see walls so spins diffuse freely; when t gets larger, 
it gets into the fast diffusion limit {ADt » (L/2)^), 
where spins see the wall frequently. As diffusion is more 
restricted in the fast diffusion limit, the effective diffu- 
sion speed is smaller than that in the free diffusion. As 
a result, the relaxation rate in the fast diffusion limit is 
smaller, similar to the explanation of the edge enhance- 
ment effect. 

We numerically evaluate 5R[/r(ajo)] as a function of 
time to reveal how l/Tp changes from the slow diffu- 
sion limit to the fast diffusion limit (Fig. 0]). Values 
of the parameters used in the evaluation are assigned as 
(x') = 0, D = 1 cmVs, L = 1 cm and Wo = 1 to 1000 
rad/s. In the figure, the quantity loq^[Jj;{u}o)] is actu- 
ally plotted for the purpose of comparison. Therefore, in 
the slow diffusion limit, ujq^[Jx{ujo)] = 2D; whereas, in 
the fast diffusion high pressure limit, ujI^[Jj;{ujo)] = D. 
As shown in the figure, when t is small, i.e. in the slow 
diffusion limit, the relaxation rate oscillates around 2D. 
When t becomes larger, the oscillating amplitude of the 
relaxation rate decreases and the mean of the oscilla- 
tion converges to the fast diffusion results. The final 
value of the fast diffusion result depends on the ratio of 
Td/ri, see Sec. II B. When uiq — 1000, it is in the high 
pressure limit and uj^JxiiOo) converges to D, which is 1 
cm^/s in our case; and when uig = 1 rad/s, it is in the 
low pressure limit and ujqJx{ujo) = ujqL'^/120D = 1/120 
cm^/s, see Eq. (p9|) . The characteristic time to distin- 
guish the slow diffusion limit from the fast diffusion limit 
is also Tc = (L/2)2/4L) = 0.0625 s. As Tc is usually 
small in practice, Tps measured by experiments are usu- 
ally in the fast diffusion limit. Nevertheless, When D is 
small enough or alternatively the cell dimension is large 
enough, the characteristic time Tc can be rather large and 
it is possible to measure the longitudinal relaxation rate 
in the slow diffusion regime. 



V. CONCLUSIONS 

A new approach based on Redfield theory is developed 
to calculate magnetic field gradient-induced longitudinal 
and transverse relaxations of '^He gas. As an extension 
to the method developed by McGregor, the newly devel- 
oped approach works in all diffusion regimes, including 
the intermediate regime. It can also explain the edge 
enhancement effect in the slow and intermediate diffu- 
sion regime, which shows the ability to capture all the 
relaxation related behaviors in one single unified model. 
It also has an advantage in terms of numerical Simula- 
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FIG. 4: (Color online)The curve with ujq = 1000 rad/s is 
definitely in the high pressure limit, the normalized relaxation 
rate evolves into D {D = 1 cm^/sec) as expected. The inset 
figure shows the relaxation rate when cjq = 1 rad/s, which is 
in the low pressure limit. It evolves into loqL* /120D, which 
is 1/120 cm^/s, when t becomes large. The other two curves 
are in the intermediate region. 



tions, because one only needs to re-compute the prob- 
ability density function for new geometries. Since the 
density function is easy to obtain, the relaxation rates 
can be computed readily as described in this manuscript. 
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Appendix A: Simplification of the Position 
Autocorrelation Function in the Slow Diffusion Limit 



In the slow diffusion hmit, iDt « (L/2)2, Eq. ^ 
is simplified by expanding all exponential terms in terms 
of t, to and t to the first order. Trigonometric functions, 

such as sin( — i—^), can be expanded in terms of , 



the first order. The first term of Eq. (|13p becomes 






n=l,3 



^ = Dt. 



(Al) 



becomes 

4 {x') L 
7r2 

4(a;' 



El .mr.mr(x'), n^ir'^Dt. 
-ism( — )^^(1 y^-) 



n=l,3 

, 2 oo 
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2 ' L '' L2 
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The fourth term is the same as the third term, which 
is also {x') ; and the last term becomes 

16L2 ^ ^ /fc2 



E E 



n 
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,k'K , 
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x(l- 



t?t:'^Dt pTT^Dto kV {x'f 



L2 



2i2 



). (A3) 



The first term in the bracket of Eq. (|A3p is evaluated to 
converge to — L2/12. The evaluation of the second term 
yields Dt. The third term and fourth term are the same, 
except for different prefactors. They are evaluated to be 
2_Dio and {x') , respectively. Collecting all these terms 
together, the autocorrelation function of x becomes 



ix{t-T)-{x'))ix{t)-{x')) 
_L2 

~ 12 



Dt +{x'f -{x'f -{x'f 



12 



1st term 
2 



2nd term 3rd and Ath term 



+ Dt + 2Dtn + {x'y +0{1/L) + 0{f) 



bth term 

2D{t~T). 



(A4) 



An alternative way to obtain Eq. (|A4[) is to solve the 
diffusion equation in free space as slow diffusion is equiv- 
alent to free diffusion. In this case, the conditional prob- 
ability function p{x,t\xo,to) is known to be 



p{x,t\xo,to) 



1 



y/4:TTDi 



(A5) 



In the free diffusion, the diffusion equation as well as 
the autocorrelation is translational invariant. Therefore, 



ill to (x(i - r) - {x')){x{t) - {x')) = x{t - T)x{t) and 



x{t — T)x{t) 



dx\ 



xoe 



i^-^or 



-dx 



y/iirDto J-oo VAttDt 
= 2D{t - r), (A6) 



The second term {x') is unchanged, and the third term which is exactly the same as Eq. (|A4|) . 
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